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Ordinary muon capture on a proton in manifestly Lorentz invariant baryon chiral 

perturbation theory 

Shung-ichi Ando*^ and Harold W. Fearing* 
TRIUMF, Vancouver, British Columbia, V6T 2 A3, Canada 

The amplitude for ordinary muon capture on the proton is evaluated, through the first four orders 
in the expansion parameter, in a manifestly Lorentz invariant form of baryon chiral perturbation 
theory. Expressions for the low energy constants in terms of physical quantities are obtained in 
each of the several renormalization schemes which have been proposed for forcing the relativistic 
approach to obey the same counting rules as obtained in heavy baryon chiral perturbation theory. 
The advantages and disadvantages of these schemes are discussed, using the muon capture results as 
an example, with the aim of gaining insight as to which scheme is preferable for practical calculations. 

PACS numbers: 12.39.Fe, 11.30.Rd, 23.40.-s, 13.60.-r 

I. INTRODUCTION 

<: 

[~~. ' Ordinary, or non radiative, muon capture (OMC) has always been an interesting process because, unlike beta decay, 
7-H , there is a sufficiently high momentum transfer to explore weak-nucleon form factors away from q 2 = 0, where g M is the 
four momentum transfer between lepton and hadron currents. The induced pseudoscalar form factor, Gp, has been of 
particular interest and OMC is the main source of information on this form factor, see, e. g., Refs. Originally 
the amplitude for OMC was written down as the most general form of the current and OMC was simply an empirical 
way to determine the coefficients of this most general current. More recently it has become possible to calculate 
these form factors from a more fundamental point of view using an effective Lagrangian of chiral perturbation theory 
(ChPT) which incorporates the symmetries of QCD Such calculations for OMC HJISGJ or f° r the electromagnetic 
form factors of the nucleoli have been carried out in so called heavy baryon ChPT (HBChPT) which involves a 
Foldy-Wouthuysen-like expansion of the Lagrangian in powers of the inverse nucleon mass. A comprehensive, modern 
review with references to earlier work can be found in 

The possibility of a fully relativistic ChPT approach to OMC, or any other process, has been elusive until recently. 
Qh! However there now have been several suggestions 0, [n|, ^2 EJ El for relativistic approaches 0] . The difficulty 
with relativistic theories has been the fact, as pointed out by Gasser, et al. 16], that they do not obey a counting 
' procedure which would allow one to associate multiple loop diagrams with higher powers in an expansion in a small 
parameter. In fact in the relativistic approach these multiple loop diagrams, specifically those involving the relativistic 
nucleon propagator, can contribute to lower orders and so there is not a well defined prescription for deciding what 
diagrams to keep. One solution to this problem was the HBChPT approach in which diagrams involving more and 
more loops contribute at higher and higher orders in an expansion of the amplitude in powers of a typical momentum 
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scale divided by the nucleon mass. While this approach works, it has the disadvantage of not being manifestly Lorentz 
invariant and of requiring increasingly complicated vertices as the order increases. The recent proposals for relativistic 
ChPT resolve this problem in a different way by showing that it is possible to define renormalization procedures for a 
manifestly Lorentz invariant theory which generate the same type of counting scheme which is present in HBChPT. 

This paper thus has a number of aims. OMC is one of the simplest non trivial processes, and so we want to use 
it as laboratory to understand how the relativistic approach and the various renormalization schemes are applied to 
a practical case. We also want to compare and contrast the various proposed schemes to see if one is preferable for 
detailed calculations or if there are alternative methods which achieve the same result but which are easier to use. 

We also want to obtain the muon capture amplitude and expressions for the low energy constants (LEC's) in a 
unified and consistent approach. In the relativistic approach one can obtain these quantities to one higher order than is 
easily possible in HBChPT. This results from the fact that the Lagrangian, and the vertex operators originating from 
it, increases in complexity with increasing order much more rapidly for HBChPT than for the relativistic approach. 
The OMC amplitude accesses the weak nucleon currents, vector and axial vector, the former being essentially the 
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same as the electromagnetic current. There is enough information available to determine all of the LEC's which 
appear. However, just as in HBChPT, the only unused information, once the LEC's are determined, serves only to 
give the well known expression for Gp in terms of the pion-nucleon coupling and the axial radius. Thus this evaluation 
of the OMC amplitude is mainly a way of determining the LEC's and does not provide a 'new' number for the rate. 

Although many of the pieces have been obtained before in separate calculations, |16l ll 7L list l2Cl l2l| the results 
for the LEC's depend on the details of the calculation and it becomes dangerous to lift these values for the LEC's from 
disparate calculations and use them for other processes. Therefore it is important to have a consistent, consolidated 
and practical approach as this will provide the basis for determining the LEC's which will be necessary for future 
calculations. Thus we see this calculation as a basis or starting point for planned similar consistent applications of this 
approach to processes such as irp — > rvy where there are puzzles in the HBChPT approach arising from the appearance 
of unnaturally large LEC's |2^. It is also intended to be a starting point for a similar calculation of radiative muon 
ca pture, u + p — > n + v + (2^, [24|, where there are still unresolved problems relative to the extraction of Gp 

II. WEAK FORM FACTORS OF THE NUCLEON CURRENT 

The S-matrix amplitude for the OMC process, fi+p —* n + is, with momenta defined by pn +Pi = pf + p v , is given 
in the notation of [23 by 

M = ~ lG ^Y ud u{p„) la (l - j 5 )u(p^)u(p f )T„ [V a - A a ] u{ Pl ) (1) 
where the vector and axial vector nucleon current operators are given by 

A a = G A (g 2 ) 7 Q 75 + ^V 75 . (2) 

Here Gp is the Fermi constant as obtained from muon decay, V u d is an element of the CKM matrix, m„ is the physical 
muon mass, m^r is the average of physical neutron and proton masses, — (mr , + m„)/2, and r_ is the isospin 
lowering operator, < n|r_|p >= 1. Here we do not include radiative corrections |29j and have neglected possible 
second class currents. 

Gy(q 2 ), Gm(<Z 2 ), Ga(<? 2 ) and Gp(q 2 ) are the weak form factors of interest, with the four momentum transfer 
q^ L =Pf—Pi- The vector and weak magnetism form factors, Gy and Gm respectively, are related to the isovector 
electromagnetic form factors of the nucleon. The axial form factor at q 2 — 0, Ga(0), is most accurately determined 
from neutron beta decay and Gp is the induced pseudoscalar form factor which is accurately predicted by chiral 
symmetry. All of these form factors are functions of the four momentum transfer q 2 which for OMC on the proton is 
given by 

2 -m^m 2 -ml + m^mp) 

q -> £ ■ = -0.88m (3) 

m p + m^ 

Note that we have normalized these form factors using the physical masses m M and tun and that we have used the 
ChPT sign convention for Ga and Gp, which makes them positive, in contrast to the convention which has been used 
historically and which is still used in the Particle Data Group listings [3(j ■ 



III. EFFECTIVE CHIRAL LAGRANGIAN 



A. Effective Lagrangian 



In the usual ChPT approach the effective Lagrangian is expanded in powers of a typical momentum - for this 
problem the muon mass m^, the pion mass m w , or the four momentum transfer squared q 2 - divided by a typical 
hadronic scale which we take as the physical nucleon mass. At each order the most general Lagrangian satisfying the 
symmetries of QCD is determined. 
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For this calculation we work in SU(2)xSU(2) and use for the chiral Lagrangian C x = C^n + £tt where C^n and 
C v are respectively the Lagrangians in the pion-nucleon and pion sectors. The pion-nucleon Lagrangian is expanded 
in terms of small quantities, 

C nN = £* N + 4jv + C* N + C^ N H , (4) 

where the ellipsis represents the higher order terms and the superscript denotes the order of the Lagrangian. 
The lowest order Lagrangian is given by the standard form 

4it = *(* P - m + \g A il^ • (5) 
Here the pion and nucleon fields are collected as 

tf = U = t #*/* *=(*)• (6) 

The covariant derivative D M , when acting on things transforming as nucleon fields, is defined as 

D^=(d fl + T^, (7) 

with 

= \[u\d^ - ir^)u + u{d^ - it^xJ] (8) 

and with and r M constructed from the external vector and axial vector currents as £^ — — and = + a M . 
We also have 

Up = £[iJ(£?p - ir^)u - u(d^ - U^u'] . (9) 

The parameters appearing in this lowest order Lagrangian, m,Fo,gA are respectively the 'bare' or unrenormalized 
values of the nucleon mass, the pion decay coupling and Ga(0) and the fields are the 'bare' fields. 
The higher order Lagrangians are given by 

7 23 118 

4*=I>^i*> ^^w,^ 4^ = 5^*0,*, (io) 

z—l i—1 i—1 

where the Cj, di, are the LEC's and where the basis functions Oi are given in Ref. 31] Tables III, IV and V. Note 
however that we always normalize the LEC's to the physical mass tun rather than to m. 

We will also need the Lagrangian in the purely meson sector, as the results depend on that choice as well. We will 
take the standard one 

z: w = 4 2) + 4 4) + ---, (ii) 

where the lowest order Lagrangian is given by 

4 2) = ^f< D^UiD^Uy >+^f< xU f + U X f >- (12) 
Here the covariant derivative acting on quantities transforming as U is given by 

D fi U = d ll U-ir fl U + iU£„. (13) 

We also have 

X = 2B Q (s + ip), (14) 

where s and p are the external scalar and pseudoscalar current and the chiral symmetry breaking is introduced as 
usual by taking p — and s — m with rh = (m u +md)/2 the average of up and down quark masses. The parameter Bq 
is given in terms of the lowest order pion mass (cf. Eq. IjBljl - ) by 2B$m — m?L. For the fourth order Lagrangian 
we take the Gasser-Sainio-Svarc [l6| form of the Gasser-Leutwyler Lagrangian |32j given explicitly for example by 
Eq. (D.13) of Ref. 0. 
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B. Counting rules and 'non counting' terms from loop diagrams 

In the usual HBChPT approach the Lagrangians £^ r ,£^ r ,£^ r ,£^ r , or more precisely the HBChPT expansions 
of these Lagrangians, contribute tree level diagrams respectively of order p,p 2 ,p 3 , p 4 where we mean by p the generic 
small expansion parameter, e.g., m n /mpj. As a consequence of using the dimensional regularization procedure for 

regularizing the integrals, |3^| the one loop graphs containing only vertices from contribute at 0(p 3 ) and those 

containing one vertex from £^jv contribute at 0(p 4 ). Two or more loop graphs contribute only at 0(p 5 ) or higher. 

In the relativistic approach however the counting breaks down, with multi-loop graphs contributing to lower orders, 
0(p 2 ) 7 0(p 3 ), etc. [16|, than that obtained in HBChPT. Thus one needs to develop some different scheme for ordering 
the various contributions and determining which to keep. 

From a purely practical point of view it is rarely possible, or necessary, to consider more than one loop. Furthermore 
one of the most important general results of Becher and Leutwyler 12] or more particularly Fuchs, et al. [lfj is that 
it is always possible to absorb into the LEC's those terms from multi-loop diagrams which do not obey the HBChPT 
counting rules. Thus we will here consider only one loop diagrams and assume that all contributions from multi-loop 
diagrams arc cither 0(p 5 ) or, in accord with the general result, have been absorbed in the LEC's in the original 
Lagrangians as defined in Eqs. © and i JTUj l. 

We note also, as will be seen from the explicit expressions below, that associated with each loop is the factor 
l/(47ri 7 b) 2 . By neglecting multi-loop diagrams we are neglecting terms containing higher powers of this factor, e.g. 
l/(47ri 7 'o) 4 . Such higher powers will also arise from the expansion of unrenormalized quantities in one loop diagrams 
about renormalized values. We will always drop such higher powers of l/(47ri<o) , arguing that such approximation 
is consistent with our neglect of multi-loop diagrams. 

Thus to summarize, we will work consistently to one loop, and to 0(p A ). This means that we will keep all terms 
of 0(p ) or lower except for those originating in multi-loop diagrams, which we assume to have been absorbed in the 
LEC's. We also drop terms involving 1/(4^^) to the fourth or higher power, an approximation which is consistent 
with the neglect of multi-loop contributions, which have these same factors. Note also that those one loop diagrams 
which have higher order Lagrangians at the vertices and thus which are of 0(p 5 ) or higher in the HBChPT sense will 
also be dropped, again assuming that any lower order terms from these diagrams which do not obey this counting 
have been absorbed in the LEC's. 

This choice of diagrams and terms to keep is to some extent dictated by the practicalities of doing such calculations. 
Two loop diagrams and those with many higher order vertices are difficult to handle, and normally would be considered 
only if there were some special circumstance which suggests that they would be large. Perhaps the most important 
result arising from the work of Refs. 0,^2 * s that this is a consistent procedure, i.e. that lower order contributions 
which originate in these higher order diagrams which we must neglect can in fact be absorbed in the LEC's in a way 
that preserves the usual HBChPT counting procedures. 

In the relativistic approach there will also be terms, a finite number of them, coming from the one loop diagrams 
we keep, which do not obey HBChPT counting, e.g. those of 0(p 2 ). We will flag these terms, but for now keep all 
of them explicitly, until we discuss the various renormalization schemes which have been proposed, as these schemes 
differ primarily in how they treat these 'non counting' terms. 

IV. EVALUATION OF THE NN 3-POINT VERTEX : VECTOR AND AXIAL VECTOR CURRENTS 

A. Preliminaries 

We now proceed to evaluate the amplitude of Eq. Q). The amplitude originates in the usual current-current 
interaction which couples the lepton current to the weak nucleon current. The lepton current is given by the relativistic 
tree-level current, l a = u u j a (l — 75)^^. The (V — A) weak nucleon current is calculated from the effective Lagrangian 
using the approach described above. 

There are three contributions to the weak-nucleon-nucleon vertex. The first two of these involve coupling of the 
nucleons to an external vector field and to an external axial field. The diagrams which contribute to these two 
contributions are given in Fig.^ The third contribution, to be discussed in the next section, involves coupling to the 

pion, which by virtue of its coupling to the leptonic current contributes to the overall axial weak nucleon current. 

(s) -* -* 

We take for the external vector current Vu — > + t • Vu, i.e., we divide the current into isoscalar and isovector 
part. Only the isovector part contributes to the weak current, but we will keep both for completeness, and to allow 
evaluation of some of the LEC's via a connection to the isoscalar electromagnetic form factors of the nucleon. Similarly 

(s) . . . (s) 

we take a M — > a/, + f ■ a M , but in this case will drop the isoscalar axial current a/, . 
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B. Tree level diagrams 

The tree level contributions to the amplitude correspond to diagrams 1-4 in Fig. ^ and are given by 



Miv = tyZ N * f (vp +f-v, l ) 1 ^ l ^Z N , (15) 

M 1A = igAv^VfT-a^j^tyfz^, (16) 

M 2V = i^Z^y /-^-^[(c 6 + 2c 7 )v^ +c 6 f -v^iy/z^, (17) 

M 2A = 0, (18) 

M 3V = WftfcT - gV)(p< +PM— vi s) + —T-v^i , (19) 

M 3A = mf^^[Aml^d 1 eT-a^ + d22{q 2 9^-%qu)r-a v ]^ i , (20) 

M 4y = il' / icr^q^[4(g 2 e54-4m^ 7r eio5)^ s) +2( 9 2 e 74 -4m^ei 06 )f-i/ AI ]* l , (21) 

M iA = 0. (22) 



Here the subscript V or A refers to coupling to vector or axial current respectively and the number refers to the 
particular diagram in Fig.^ For these and subsequent amplitude expressions % j and \&i are to be interpreted as wave 
functions, rather than the fields of the original Lagrangian. That is, they are still two component objects in isospin 
space, but made up of spinors u(pf) and it(pj) rather than fields. In the HBChPT counting system these tree level 
diagrams contribute to 0(p), 0(p 2 ), 0(p 3 ),0(p 4 ) respectively. Note that the nucleon wave function renormalization 
factor Zjv appears only in the Mi and M 2 amplitudes which is a consequence of the fact, as shown in Eq. i|B21(l of 
Appendix iBl that the leading corrections to Zn are two orders higher, Z^ = 1 + 0(p 2 ). 



C. Leading one loop diagrams 

The next set of diagrams consists of those one loop diagrams with all vertices coming from £^y. In the HBChPT 
sense these all contribute first at 0(p 3 ) , but in the relativistic approach they will also contribute some 'non counting' 
terms of 0(p 2 ) as well as relativistic corrections of 0(p 4 ) and higher. 

We express these amplitudes in terms of a general loop integral Iit-k....nn.... which is defined in detail in Appendix 
lAl For present purposes it is sufficient to note that Iinr....NN....[ki, ■■■■Pj, A] refers to a loop integral which contains 
a pion propagator for each subscript 7r with momenta of the form k% + 1 and a nucleon propagator denominator for 
each subscript N with momenta of the form pj + I. The loop integration variable £ [34( is chosen so that the first 
pion momentum k\ is zero and that argument is normally not put in explicitly. A is whatever is in the numerator. 

In terms of these loop integrals the lowest order one loop amplitudes are given by 

■ 2 

M *v = j£*J&»V> -7-^)I«NN\puPf, h^f+ m)^{^+ ?+ m) , (23) 

. 3 

M ^ = -T^r-^NN[pi,PfJl^f+ $+m)^ l5 (p' i + ^+m) fhfs^i, (24) 



+WK,7 A '75(# i + ?+m) frs]}*i, (25) 

M SA = ^VfT-a^NlPf, h^f+ ^+m)^]+I vN \p h ^(M ?+m) #75]}*,, (26) 

Muv = —^*ff-v^Ml], (27) 

M 11A = - % M#fi t -3 li <f<y 6 *iI v [l], (28) 

• 2 

M 13V = --U^fr- v^ N [-q, Pi ,(2e-qr^- fl-ftiM ?+m) (29) 

M 13A = 0, (30) 



6 



M ^v = izL^fr-v^l-q, {21- g)"(2 f- rf)]*, , (31) 
M 15A = 0. (32) 

D. Further one loop diagrams 

The next class of one loop diagrams consists of those with one vertex from Cf^fL with all others from C^L . In the 
HBChPT sense these would be of C(p 4 ), but again in the relativistic formulation they will have 'non counting' terms 
of lower order. 

These amplitudes are given by 

■ 2 

M *v = ^%[3(c6+2 C7 >W-c 6 f-iv] 

^ 

xI«NN\pi,Pf, Mtff+ f+ m)—-^(&+ f+ m) , (33) 

Zitln 

M 6A = 0, (34) 
M gv = 0, (35) 

M 9A = - l -B^fT-aAlMPfJl^ S + ?+m)[^-(p? Pi -l+(p f + ir(p f +£)■£) 

+ 8C3^ - 2ia^(4c 4 £ u + —q u )}} 

m N 

+ I*N\pi, [^-(pfa -£+{pi + iy i {p l + £)■£) + 8c 3 F 
m N 

+ 2ia^(4cdu ~ —Qy)]{M ?+ m) fj 5 ]\ , (36) 
m N J 

M w = J^t*M s) +f-v ll )(p i +p f ) v I^t]^ i + -J*^ s o^q v T- v^M\\ , (37) 
r m N Z£ m,N 

M l2A = 0, (38) 
2ca — 

M 16V = — E i^fa ,M ' / q v T-v a I^[-q,{2£-q) a e^i, (39) 

M 16A = 0. (40) 

E. Mass insertion terms 

The final set of amplitudes arises from the mass insertions on internal nucleon lines. These insertions come from 
the NN two point term in £^y, namely \P fAciiriQ^i. The relevant amplitudes, those for diagrams 7,10,14 of Fig. Q 
can be obtained from the underlying diagrams 5,8,13 respectively by the substitution for each nucleon propagator in 
turn 

. * , , (4i Cl mgJ . * , , . (41) 



p — m + it p — m 

An alternative procedure is to observe that m appears in the loop integrals only in the propagators, since we 
normalized all the constants in the Lagrangian to nijv not m. Thus we can use the fact that 

i(4zcimn 7r ) T ^-(- r ) = -; -(4dcimL)-. r- (42) 

This allows us to obtain the amplitudes with mass insertions by taking derivatives of the corresponding amplitudes 
without insertions. As detailed in Appendix ^ this approach, while exact, may not be as useful as one might expect 
because the derivative in effect reduces the power of small expansion parameter. Thus in some cases one has to expand 
the initial integral to higher order than needed for it alone so as to get the mass insertion diagrams to the appropriate 
order. We actually calculated these mass insertion diagrams explicitly and used this derivative procedure to check 
the results. 
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Another approach is to observe that (see Appendix [BJi the 'bare' nucleon mass which appears in the original 
Lagrangian is related to the physical mass by the relation 

m N = m-4:C 1 nil 7r + 0(p i ). (43) 

Thus a propagator with mass m —* can be expanded as 

-, ; > -, - + -, -(Ucxml n )-. - + .... (44) 

p — niN + le p — m + ie p — m + ie p — m + ie 

Thus if we replace m — » m^r in the propagators in the one loop diagrams we effectively include the mass insertion 
diagrams 7,10,14 of Fig. 2] At the same time we reduce the number of separate diagrams to be calculated and reduce 
by one the maximum number of propagators involved in the loop integrals which must be calculated. Both of these 
offer significant calculational advantages. However this expansion only works to first order, so one must keep only 
terms in the expansion of the propagators which are linear in c\ . This requires extreme care since c\ appears also in 
other places in the amplitudes so that in general there are legitimate terms involving c\ which must be kept. 

Note that simply replacing m — > mjy in all the propagators is not exactly equivalent to the direct calculation of the 
mass insertion diagrams. The expansion of Eq. Q44|l only works to lowest order, so with such a substitution there will 
be some spurious higher order terms implicitly included. Also implicitly included will be diagrams like those of Fig. [21 

(2) 

which involve two mass insertions or which involve one mass insertion plus a vertex from C]-N- These would not have 
been included in a direct calculation of the mass insertion diagrams as they would have been nominally of too high 
order. To the order we are considering most of these extra terms can be neglected. In fact in HBChPT all would be 
of higher order. However in the relativistic approach there will be a few terms arising from the 'non counting' terms 
from diagrams such as those of Fig. [2] which will appear in the amplitude in this approximation and not in the explicit 
calculation. 

To repeat, our approach was to calculate the mass insertion diagrams explicitly and thus our results may differ 
from calculations which use one of the above approximations. 

F. Summary 

By summing all of the amplitudes given above we obtain the complete weak nucleon-nucleon amplitude arising from 
interaction with external vector and axial vector fields. 

V. EVALUATION OF THE NN 3-POINT VERTEX : PION TERMS 

The third contribution which must be evaluated comes from the pion-nucleon-nucleon vertex and will contribute 
to the axial current. The diagrams which are needed are given in Fig. [3] and the amplitudes associated with those 
diagrams are given by 

Mi* = --j^y/Z^fT-T? faViy/Z^y/z^, (45) 
M 27r = 0, (46) 

2 

M 37r = ^(di8-2d 16 )%f-7?# 7 g^, (47) 

^0 

M 47T = 0, (48) 
Q 3 — 

M ^ = ^3^f^-^NN\pi,Pf, M0f+ f+m) tfy 5 (jfc+ h^i i (49) 

Me, = 0, (50) 

+ U N \pi,U+ ?+m) /7s]}*,, (51) 

-2^3*/^ ' nUtNlPf, #7s(^/+ f + m)[-4cimj^ 

-- t-((p/ + e ) ■ £ (pf + e)-o.+Pi-epi-q)- ^ 3 e ■ q - 2^4^%^]] 

m N 



Mo 
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-I*n\Pu [-4cim^ + -T-((pi + £) ■ *(p< +e)-q + Pf-epf-q) + 2c 3 £ ■ q 
m N 

~2i Ci a^e^q u }(M ?+m) ? l5 ]}y t , (52) 
M Uw = ^3*/r-ff ^75*A[1], (53) 
M 12w - 0. (54) 
Again the mass insertion diagrams 7 and 10 are obtained by making the replacement of Eq. I|41l) in diagrams 5 and 

8. 

To get the contribution to the weak nucleon-nucleon axial current from this ttNN amplitude we make the replace- 
ment (for q 2 ^ m 2 ) f • 7? — > 2iF 7T q fJ "r- a^j (q 2 — m 2 ). This arises from the addition of a pion propagator and pion decay 
vertex to the irNN vertex. Note that the parameters and are the physical ones. Since we have associated 
a V 'Zn with the amplitude Mi OT so that it is renormalized, we need to use here the renormalized propagator and 
renormalized pion decay vertex to make the overall amplitude renormalized. 



VI. EVALUATION, REGULARIZATION AND RENORMALIZATION 

The first step in evaluating these amplitudes is to reduce the numerators of the loop integrals. This is done 
using the standard algebra of Dirac matrices and the usual tensor decomposition of integrals with 1^,1^1",,.. 
in the numerator. The end result is that the full amplitude can be expressed in terms of the following loop 
integrals with unit numerator: I^[l] 1 I^ N [p,l]J NN [pi,p } ,l]J^[q,l]J^ NN [p i ,p } ,l]J^ N [~q 1 p,l], where p can 
be Pi or pf. For the diagrams with the mass insertion put in explicitly we need the additional integrals 
InNN\Pi,PiA],lNN\Pi,PiA],lNNN\PhPi,PfA],I™NN[-q,Pi,PiA] plus the corresponding ones with the roles of p { 
and pj interchanged. 

In HBChPT these integrals are evaluated using dimensional regularization to extract the divergences, which are 
then absorbed in the renormalization of the LEC's. In the relativistic approach this procedure works in the same way 
for integrals involving only pion propagators. Thus for example we have in standard fashion for dimension d ~ 4 

iAn = ^[R + H^)} (55) 

where 

1 4 — d 

R = +7-l-ln(47r), e = — ^— (56) 

with 7 = -r'(l) = 0.577 

In the relativistic approach the same procedure applied to integrals containing nucleon propagators leads to the 
MS scheme in which the R's are all absorbed in the LEC's. The amplitudes however still contain finite terms which 
do not obey the usual counting rules of HBChPT. Thus for example one loop diagrams which are nominally 0(p 3 ) 
may contain contributions at 0(p 2 ) and likewise those nominally 0(p 4 ) may contain also terms of 0(p 2 ) and 0(p 3 ). 
There have been two somewhat different, but similar, methods proposed to resolve this problem. In the 'Infrared 
Renormalization' (IR) scheme proposed by Becher and Leutwyler |12| the loop integrals are divided into two parts. 
An 'infrared singular' part contains non integer powers of the small expansion parameter and a 'regular' part contains 
only integer powers. They then renormalize the integrals by dropping the regular part and absorbing the infinities 
of the singular part, i.e. the terms proportional to R, in a renormalization of the LEC's. Note that 'drop' means 
'absorb in the LEC's' or equivalently 'cancel via counterterms in the Lagrangian'. Thus in this approach the infinities 
which appear in both singular and regular parts are in effect combined and absorbed in the LEC's in the same way 
as would be done in the MS scheme. The difference arises in that in the IR scheme additional regular polynomial 
terms, including all those which do not satisfy usual HBChPT counting, are also absorbed in the LEC's. 

The other approach, the Extended On Mass Shell (EOMS) scheme of the Mainz group [T^j, first uses the usual 
dimensional regularization to extract the terms proportional to i?, which are then absorbed in renormalizations of 
the LEC's in exactly the same fashion as in the MS scheme. In the second step the amplitude for each individual 
diagram is examined and those terms, all polynomials in the expansion parameter, which do not obey the counting 
rules, as used in HBChPT, are determined. This finite set of terms, the 'non-counting' terms, are then dropped, i.e. 
absorbed into the LEC's. This approach thus eliminates a somewhat smaller set of terms than does the IR approach. 

In practical applications the EOMS approach involves a number of subtleties. These subtleties, basically amounting 
to choices of conventions, affect the specific terms absorbed in the LEC's and thus make little difference as long as 
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one is considering just one process. They simply change slightly the numerical values of the LEC's. However if one 
wants a consistent scheme to be applied to a variety of processes, as is our intent here, it is necessary to discuss the 
various choices and to define exactly what conventions we take, as one must use the same conventions in subsequent 
calculations or when using values of the LEC's extracted by others. 

First, when extracting the non-counting terms from each diagram it makes a difference whether one first expresses 

the amplitudes as a function of the original mass m appearing in the Lagrangian, the mass in the chiral limit m, or the 
physical mass m^. Differences are 0{m 2 .) and would thus be higher order corrections in the HBChPT scheme. In the 
relativistic approach however such corrections to, say, 0(p 2 ) non-counting terms can enter as 0(p ) terms which are 
kept. Thus absorbing a non-counting term expressed as a function of m leads to a slightly different renormalization 
than would be obtained by absorbing the equivalent term expressed as a function of ■ Here we always express the 
amplitudes in terms of the physical mass tun before isolating the non-counting terms. 

Similarly the log terms in the amplitudes can be expressed as ln(m 2 //z 2 ) which makes /i = m the logical choice, 
since then these terms vanish, or as \u(m 2 N / [i 2 ) which makes /i = tojv the logical choice. We have kept the scale 
parameter fi explicit until the end, but have used uin in the amplitudes and eventually for fi. 

An essentially similar effect arises in the approach used for including the diagrams with mass insertions on internal 
nucleon lines. We evaluated each diagram explicitly so that for each of the 0(p 3 ) diagrams with an internal nucleon 
line, i. e., Fig. ^diagrams 5,8,13 and Fig. |21 diagrams 5,8, there is an associated diagram of 0(p 4 ), Fig. ^diagrams 
7,10,14 and Fig. [31 diagrams 7,10 respectively. We then looked at each diagram in the two sets to determine the 
non-counting terms for each. An alternative approach used in |18| replaces m — > in the 0(p 3 ) diagrams and then 
later expands to first order in the difference to get the mass insertion contributions. In this approach a term of 0(p 2 ) 
before expansion would be dropped, as the original diagrams are nominally of 0(p 3 ). However had the expansion 
been done first, the expansion of such terms would give pieces of 0{p A ) which one would keep, and which in fact are 
some of those arising in the diagrams with explicit mass insertions. 

Finally observe that, although the EOMS scheme is based on extracting from each individual diagram those terms 
which do not obey the nominal counting rules and absorbing those terms in LEC's, that procedure does not ensure 
that each individual diagram obeys the counting rules. The exact same statement can be made for the IR scheme 
and a very similar statement can be made for the MS scheme, where the renormalization of the LEC's to remove the 
infinities ensures that the amplitude is finite, but not that each individual diagram is finite. In all these cases the 
situation occurs because in general there are not always LEC's available to absorb terms from individual diagrams. 
The simplest example of this can be seen in the calculation of the q 2 = limit of Gv ■ In the EOMS scheme explicit 
calculation shows that the diagrams 7,8,10,13, and 14 of Fig. ^all contribute non-counting terms. Many of these are 
removed by the renormalization provided by Zn (appearing in the tree level diagrams 1 and 2) as given in Eq. (|_B21|) 
of Appendix^] However two contributions remain, those from diagrams 8 and 13. There are no LEC's available here 
and so no way to absorb these as individual terms. Instead what happens is that these two contributions cancel each 
other so that the sum of the amplitudes from all individual diagrams contains no non-counting terms. 

Similarly in the MS scheme diagrams 5,7,8,10,11,13,14, and 15 of Fig. all contain infinite terms proportional to 
R. Again Zn renormalizes some of these away, but there are a number of terms left and no available LEC's to absorb 
them. Instead they cancel among themselves. 

Naively it is perhaps obvious that this happens as individual diagrams are not physically measurable quantities 
and thus do not necessarily satisfy physically relevant constraints such as counting or finiteness. More precisely, 
the argument that the finite number of non-counting terms, or the infinities, can be absorbed in the LEC's (or in 
counterterms) relies on the fact that the Lagrangian contains all possible counterterms allowed by the symmetries. In 
this case the relevant symmetry is current conservation which, as is well known, ensures that the weak vector coupling 
(or the isovector electromagnetic coupling F\) is not renormalized by the strong interactions. This symmetry is 
obeyed by the full amplitude, but not by individual diagrams. Thus it perhaps should not be a surprise that terms 
are generated for individual diagrams which cannot individually be absorbed in counterterms. The fact that the sum 
cancels as required however is a clear check on the correctness of the calculation. 

Thus to summarize we consider three approaches, the MS, IR and EOMS schemes. In all three one first regularizes 
the integrals using dimensional regularization and renormalizes the LEC's to remove all infinities in the usual way. 
For IR and EOMS additional sets of finite terms are extracted and absorbed in the LEC's so as to preserve counting. 
The difference between the two is simply in the explicit terms extracted. 

Note that these three approaches will not lead to any different predictions for measurable quantities. The formulas 
for such quantities in terms of the LEC's will be different, but that will be compensated by different formulas for the 
LEC's in terms of unrenormalized quantities and different numerical values for the LEC's. 
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VII. RESULTS 



We thus proceed as outlined above, i.e. we extract R from each of the remaining loop integrals using the standard 
dimensional regularization as in Appendix El The integrals are first separated into two parts according to the IR 
prescription. The parts involving R are recombined and the renormalization of the LEC's to absorb R proceeds just 
as in the usual MS scheme. Thus the original LEC's in the Lagrangian, x, where x stands for any of the LEC's, 
are eliminated in favor of their MS renormalized values x r . The finite parts are then expanded in powers of the 
small parameter and terms through 0(p A ) are kept. We flag those terms originating in the infrared regular part of 
the integral with a parameter f3 IR , as detailed in Appendix El We assume that q 2 /ml < 1, but not necessarily very 
small, but for simplicity keep only terms linear in q 2 in the final results. We also replace the original parameters of 
the Lagrangian, m,mo n> Fq with their physical values as determined in Appendix 151 Then the contribution of each 
diagram is examined and those terms which do not obey counting and which would be dropped in the EOMS scheme 
are flagged with a parameter [} EOMS . 

The full amplitude is then evaluated and put in the form of Eq. Q which gives the vector and axial vector weak 
nucleon currents, V^^A^, appropriate for muon capture and allows us to extract explicit expressions for the various 
form factors in the equation. Two further renormalizations are then performed. The first expresses x r in terms of 
the EOMS renormalized quantities x EOMS and is determined by requiring that all terms flagged by p EOMS must be 
absorbed. The second expresses x EOMS in terms of the IR renormalized LEC's x IR and removes all the remaining 
terms proportional to (3 IR . The expressions for the renormalizations seem to be unique, except for the few cases 
where only a combination of LEC's appears, as long as each renormalization involves only terms with the same power 
of m v . 

The weak form factors expressed in terms of the IR renormalized LEC's are given by the following. In these relations 
we have always used physically measurable masses, m^,m n and coupling F v , have taken the scale factor [i — > mjv, 
and have kept only terms linear in q 2 , but otherwise have kept all terms consistent with an expansion of the amplitude 
to fourth order in the expansion in the small parameter. 
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The correction terms needed to renormalize the LEC's and other parameters arc all proportional to l/F 2 . Thus 
in terms already containing this factor it doesn't matter which of the various renormalizations are used as we have 
consistently neglected terms of order 1/F£. Thus to simplify the notation in the above, and also the equations below, 
we have left off the superscript, r, IR, or EOMS, on g A and the various LEC's when they appear in terms already 
containing a l/F 2 . Note however that for numerical work we will always use the appropriate LEC wherever it appears, 
i.e., when working in the IR scheme, all LEC's will be the IR values, and similarly for the other schemes. 

Since we kept the isoscalar component of the external vector field we can also obtain the isoscalar electromagnetic 
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We find for the MS 1 renormalized LEC's expressed in terms of the LEC's of the original Lagrangian: 
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Note that the renormalizations of g A given in Eqs. iJBSJ, (|77() . and (|89|1 originate from terms that survive in the 



chiral limit and thus they renormalize the original g A appearing in the Lagrangian to g A in the chiral limit, 3a- Just 
as discussed for the mass in Appendix [5] very often a counterterm to perform this renormalization is included in the 
original Lagrangian and gA is assumed from the beginning to be gA in the chiral limit. 



VIII. NUMERICAL EVALUATION OF LEC'S 



In the preceding sections we have obtained the result for the complete amplitude for OMC as expressed in Eq. 
using the values for the couplings from Eqs. (|5?|) - l|6T)|) . This amplitude is expressed in terms of the physical masses, 
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the pion decay constant F n = 92.4 MeV, the external parameters ci, 02,03,04 and the sets of LEC's x r , x IR , or 
x EOMS , depending on the case being considered, where x stands for g^, eg, C7, de, d-j, die, dig, d 2 2, 654, 674, e^s, eio6- 
To determine these parameters we use available data from measurements of weak and electromagnetic form factors. 

(v) (v) 

For the vector current we have information on the isovector form factors and F 2 1 equivalent to Gy and 

(s) (s) 

Gji, and on the isoscalar form factors f{ and F 2 . The static values of the magnetic form factors are given by 
F [ 2 v) {Q) = n p -n n and if°(0) where the proton and neutron anomalous magnetic moments are taken as 

k p = 1.7928 and n n — —1.9130. We define the slopes of the various form factors in the usual way 

F(q 2 ) = F(0)(l + ^ <r 2 >) (101) 

where q 2 is the square of the four-vector momentum transfer and < r 2 > is the rms radius. We take the values of the 
rms radii for Fi, F2 in the isoscalar and isovector cases from Mergell, et al. [j^l and thus use < r 2 ^ >= (0.765fm) 2 , < 
r 2 2 (v) >= (0.893fm) 2 ,< r 2(s) >= (0.782fm) 2 , < r 2(s) >= (0.845fm) 2 . 

Information on the axial current comes from neutron beta decay which gives Ga(0) — 1.2695 ± 0.0029 30] and 
from antineutrino-nucleon scattering [3(| which gives the axial rms radius < r 2 A >— 0.42 ± 0.04fm 2 . We use for the 
pion-nucleon coupling constant G^NNi™ 2 ,) = 13.0 ±0.1 |37| . 

There is one remaining unused equation, Eq. (|60|l . which gives the well known expression for Gp(q 2 ) in terms of 
G-n-jviv^TO 2 ) and ^22, which can be determined from < r\ >. In principle, if Gp were well measured, this could be 
used as an alternative to one of the equations to determine the LEC's. In view of the uncertainties in the experimental 
value of Gp 0] however this is best used to predict Gp or simply as a consistency check. 

Finally we need the external parameters 04,02,03,04 which can be obtained from pion nucleon scattering. One 
should in principle evaluate these via a complete calculation of pion-nucleon scattering consistent in order and in 
its details with the calculation here. That is beyond the scope of the present paper. So for present purposes we 
will simply take the results of a tree level fit obtained by Becher and Leutwyler [l7|. namely, c\ = — 0. %ra~^,c 2 = 
2. 5771^,03 = — 4.2m^ , C4 = 2. 3m^ . These parameters appear only in higher order terms, so this approximation is 
probably sufficient. 

We have identified above nine bits of experimental data to be used to evaluate the parameters. However there are 
twelve unknown parameters. Note however that at least for the IR and EOMS schemes at leading order only certain 
combinations appear. Thus we define 

~9a = .g^ + 4m 2 rf{^, 
ci R = ci R -lQmlm N e IR 6 , 

= c I 7 R ~8mlm N (2e[ R -e I 1 R ), (102) 

with an analogous definition for the EOMS and MS schemes. For the IR and EOMS schemes, which obey counting, 
the to 2 coefficient in these definitions means that we can replace all 3,4,06,07 appearing in higher order terms with 
gA,ce,cj. Thus we eliminate all instances of die, eios, an d eio6, and so have enough input data to solve uniquely for 
the nine parameters. 

For the MS scheme however this doesn't work. Because of the non counting terms the replacement 17^,06,07 — > 
gA,Ce,cj leaves some instances of (fie, eios, and eio6 which are not of higher order. Thus we need to assign values 
to these LEC's in order to solve for the others. Since there is not enough experimental information available we will 
simply try a couple of arbitrary cases to get a feel for the sensitivity of the results to these LEC's. In particular we will 
take, for a case MS-a, die = e i05 = e i06 = 0. As an alternative we will take for case MS-b, die = e ios = eio6 = 1, 
expressed in appropriate units. This latter choice is arbitrary, but should correspond to a 'natural' size for these 
LEC's. 

To actually solve for the LEC's for say the IR case we take Eqs. H57|) - (|59|) . i|tjl|) . l|t)2[l and I|tj4fl and express all of the 
LEC's in terms of their IR forms, so that the equations are expressed purely in terms of IR quantities. We then solve 
these equations self consistently, using the experimental input given above, for all the LEC's. In particular this means 
that we solve the cubic equation for g IR and use that value in the other equations to solve for the other LEC's. To get 
the EOMS case we use Eqs. I|65f) - (I76II to replace the IR LEC's with their EOMS forms, dropping higher order terms 
as appropriate, so that the equations are given entirely in terms of EOMS quantities, and then repeat the solution 
procedure. Note that this procedure corresponds to what one would do if one were using the EOMS scheme from 
the beginning. It is not quite the same as simply using Eqs. J§2J-J7BJ| to get the EOMS LEC's from the IR results 
because of the numerical consequences of higher order terms which would be treated slightly differently in those two 
approaches. Finally one gets the MS results in analogous way, though here as noted above, for that case we have to 
choose values for die, eios, and eio6- 
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TABLE I: Results for g~A and the various LEC's in each of the renormalization schemes discussed in this work. Given for 
comparison are results from |2CB and |lq | converted to account for the normalization, to mjv vs m, which we use, and in the 
case of Refs. |l8 f .|l9 | for the different combinations of cq and £7 they use. The parameters <?a, ck, and £7 are dimensionless, 
and the di and a have respectively units of GeV -2 and GeV -3 . The cases labeled MS'-a and MS-b involve arbitrary choices 
of di6,eio5, and eio6 as described in the text. 
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IR 


0.9568 
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0.30 
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EOMS 
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6.35 


-3.26 
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-0.49 


-0.17 
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0.26 
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MS-a 


-0.6244 


2.52 


-0.49 
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0.28 


2.72 
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0.0 


0.0 
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2.29 
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-0.44 


-0.48 


2.29 


0.26 
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0.25 
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-2.49 


-0.75 


-0.54 






0.19 


1.59 









The results for the LEC's obtained as described above by consistently solving all the relations available from the 
OMC amplitude are given in Table [I] together with available results obtained by others. 

First we should comment on the comparison with previous results. There have been two previous calculations of 
the electromagnetic form factors and the corresponding LEC's in relativistic formulations, Refs. [13,^3 • While our 
results are qualitatively the same, there are differences in detail. 

Perhaps the main difference in principle is the value of gA used. In both of these previous works <?a was taken to be 
gA — Ga(0) = 1.26 which is the lowest order result of Eq. Ij59(l . Also since d\Q doesn't appear explicitly in the vector 
current it was not necessary there to distinguish between gA and cja- We however expressed everything in terms of 
gA and solved Eq. (|59|l consistently to the order of the calculation to obtain a value of gA- Since gA appears in many 
places, and in particular in the corrections to all the other LEC's, this made a difference, significant in some cases, 
in the values of the LEC's obtained. In a purely formal sense the corrections to gA, i- e. the differences between gA, 
gA, and the lowest order approximation 1.26 are all of higher order. Thus in principle the use of any of these three 
interchangeably in the formulas for the LEC's would be consistent with our other approximations. The fact that it 
makes a difference simply reflects the fact the the higher order terms are not always small, i. e. that the expansion 
doesn't always converge well. However since we have the information, via Eq. to calculate the corrections to gA, 
it seems appropriate and preferable to use that information consistently in obtaining the other LEC's. Finally note 
that some further differences arise because in our self consistent solution for gA the results are different for the IR 
and EOMS cases, because die is different for those cases. 

Additional smaller differences arise because we used the rms radii appropriate to the Dirac and Pauli form factors 
F\ , F[ 2 which were the form factors calculated directly, rather than converting to radii appropriate for the Sachs form 
factors. This affects some of higher order terms and seems to affect d e particularly. Also, as discussed above, there 
are different options for including mass insertions and for expanding to get the non counting terms and we did not 
always use the same conventions as in previous work. The value of C4 used was slightly different than the one in 
Ref. |2f|. Finally we expressed everything in terms of the physical mass mjy instead of m. Again formally these 
should be interchangeable, but numerically it made a difference in some cases. 

While we think the use of the self consistent value of §a is a definite improvement in principle over previous works, 
the other differences are really just differences in the details of the calculation. The fact that they make a numerical 
difference in the values of the LEC's just reinforces the statement made at the beginning. Namely, if one wants 
values of the LEC's which can be used in further calculations one must be sure that the same approach and the same 
conventions and approximations are made. Otherwise it is dangerous to simply lift results from one calculation to use 
in another. 

Now let's look more carefully at the results of Table [I] Note that gA,£e,cr differ for the IR and EOMS cases. 
Since the underlying parameters gA,CQ,c 7 don't change (cf. Eqs. ljfj5*|) . JBSJl, and these differences must be due 

to differences in di6, eios, e\o& as given in Eqs. (|7U|l . (7SJ), and l(7tj|) . The LEC d§ is very different for IR and EOMS 
schemes, and also varies from previous results, this is apparently because of strong cancellations among terms, which 
make it very sensitive to the small corrections. Our value of cq for the EOMS case differs significantly from that of 
Ref. 01 apparently because of the c\Cq and C1C7 terms we have (cf. Eqs. (|75|l and (|76|l ) which they have not kept. 
These terms seem to originate in the different way of including the mass insertions which we used. 

In the MS scheme the parameters gA,c&,C7, and de change fairly dramatically as compared with values obtained 
in the IR or EOMS schemes. Apparently de is still sensitive to cancellations and the other three contain large non 
counting terms, which also do not vanish in the chiral limit. Had we adopted the common procedure of first adding a 
counterterm to the Lagrangian to renormalize gA to the chiral limit, such large terms would not be there for gA, and 
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presumedly it would be the same as for the IR and EOMS schemes. However such large terms would still be present 
for cq and c-j and would still affect c/a, c&, C7 via the values of d Wl e 10 5, eio6 buried in them. 

Note that all the results for the MS scheme are dependent on the somewhat arbitrary choices made for die, eios, ei06- 
The two illustrative cases correspond to values of zero for these LEC's and values of unity in natural units. Many of 
the LEC's are similar for the two cases but a few, particularly £7, change a lot. Clearly if one wants to seriously use 
the MS scheme, it will be necessary to pin down d\Q, eios, e\oe from some other process. 

Finally we should make a few general remarks. All three of these schemes, since they differ only in how they absorb 
or do not absorb the finite non counting terms in the LEC's will give the same values for the amplitudes. One might 
hope that one scheme or another would, say, lead to all small LEC's which could be neglected. This does not seem to 
be the case and there does not seem to be any general pattern emerging when we compare the three schemes. The 
parameters gA, cq, C7 are perhaps a bit smaller in the MS scheme than in the others, indicating that the specific non 
counting terms which are kept explicit in the MS scheme but absorbed in the LEC's in the other schemes are large. 
However this does not persist for the d's or e's which are of the same size, or maybe smaller, in the IR and EOMS 
schemes as in the MS scheme. 



IX. DISCUSSION OF DIFFERENT APPROACHES 

In previous sections we have described an explicit calculation - that of the amplitude for OMC - carried out in three 
different schemes for Lorentz invariant chiral perturbation theory. In this section we want to compare and contrast 
these schemes, particularly from the point of view of how best to do a practical calculation. 

First, as a matter of principle the IR and EOMS schemes are major advances in our understanding of how to handle 
Lorentz invariant ChPT calculations. Such approaches show that in general it is possible to rewrite relativistic ChPT 
so that it obeys the same counting rules as HBChPT, which thus solves the problem with such theories raised in [l6| . 
It was also shown, particularly in 0] that, unlike HBChPT, these schemes preserved the correct analytic structure 
of the amplitudes. That feature has not been important for the OMC calculation, but can be for other processes. 

Thus we now know that in a relativistic theory the choice of number of loops and the choice of the order of the 
Lagrangian to use at each vertex can be made in a rigorous way that preserves HBChPT counting, and that low 
order contributions from higher order diagrams can all be absorbed in a consistent way in the LEC's. From the point 
of view of a practical calculation that means that the choice of diagrams and vertices can be made essentially as in 
HBChPT. 

Once that choice is made however, from a practical point of view, one has options. We have considered three 
possibilities for renormalization: IR, EOMS and MS. All three treat the infinities, i. e. the R terms in the same way. 
They differ only in which subset of the set of finite terms which do not obey counting are absorbed in the LEC's. 
Thus the LEC's will have different numerical values in the three schemes and the formulas for measurable quantities 
will look different. But all three will give the same predictions for measurable results. Once the general principles 
have been used to choose the diagrams to be considered, any one of the three schemes could be used consistently for 
practical calculations and would give equivalent results. 

We can discuss however some of the pros and cons of the three schemes, relative to practical calculations. 

Consider first the IR approach. It absorbs the largest number of terms in the LEC's and as a consequence the 
formulas tend to look simpler. However one might be hiding known physics by absorbing such terms. This approach 
is probably the simplest of the three as long as one does not need to work out the exact formulas for renormalization 
of each of the LEC's. This is because if one just 'drops' the terms which would later be absorbed one can drop a 
lot of integrals - all with only nucleon propagators - and thus reduce the number of diagrams to be calculated. If 
one calculates explicit formulas for the renormalization of the LEC's, which we have done here, though it would not 
normally be really necessary, then all diagrams have to be calculated for all three schemes. 

In contrast the MS scheme absorbs none of the finite terms. It is thus closest to the historical approach of describing 
a process by a set of Feynman diagrams. Some non counting terms will appear, but may be considered to have physical 
significance. An example of this can be seen in the classical approaches to radiative corrections to neutron beta decay 
where certain terms, which in the relativistic ChPT approach seem to originate as non counting contributions from 
diagrams too high order to keep [29j, appear explicitly in the standard Feynman diagram approach |38j . have been 
discussed individually |39|. and are considered relevant. 

The MS scheme requires more effort than the IR scheme, if explicit formulas for the renormalization are not 
required, as one must always calculate all diagrams. Since there are non counting terms still present, the grouping 
of LEC's to reduce the number of independent quantities to be fitted to experiment, as done in Eq. () 111 2(1 . will not 
necessarily work, as we saw for the present calculation. This is a serious disadvantage for a single calculation as it 
increases the number of LEC's to be evaluated from data. It might be less of a problem for a series of calculations as 
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it is unlikely that the same grouping will work for all processes and so in that case for all the schemes one probably 
has to evaluate all LEC's individually anyway. 

The EOMS scheme is somewhere in between the other two. It absorbs the minimum number of terms necessary 
to get counting. It thus may preserve some of the good things about the MS scheme while still solving the counting 
problem. It however requires the most work of all as every diagram must be evaluated and then one must look at 
each diagram individually to determine which terms to subtract. It also requires a careful statement of conventions, 
as discussed above. 

In a general sense the LEC's absorb our ignorance, so it would seem that one would want to leave explicit as much 
known physics as possible, and absorb as little as possible into the LEC's. Ideally the LEC's representing unknown 
physics would then get small. This is the general philosophy behind attempts to include explicitly additional degrees 
of freedom, such as the A |4fll4ll | or vector mesons [13,E3- Thus smallness of the LEC's might be a criterion for the 
choice of scheme. One has no knowledge of the size or sign of the sum of terms contributing to an LEC from higher 
order diagrams however. Also, in the present example, OMC, there is no obvious choice leading to small LEC's, so it 
is not clear how to implement this criterion. 



X. CONCLUSIONS 



We can thus summarize as follows. We have evaluated the OMC amplitude through 0(p 4 ) in the three schemes, 
MS, EOMS and IR. Using available data we have solved self consistently for the nine LEC's which appear in the IR or 
EOMS schemes. The MS scheme requires three additional LEC's, for which further data would be required. Similar 
evaluations of the LEC's for the vector current have been done before, and our results differ from these primarily 
because we have self consistently solved the equations coming from the axial current for qa and have used that value, 
rather than the lowest order result used in previous work. Many subtleties and details of the calculation also affect 
the numerical values of the LEC's, which indicates that before using these or other values of the LEC's to calculate 
new processes it will always be necessary to make sure that the new calculation is done in exactly the same way as 
that used to extract the LEC's. 
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APPENDIX A: LOOP INTEGRALS 



We will define the general loop integral in d dimensions containing i pion propagators and j nucleon propagators 
and corresponding to the momenta as in Fig. 0] as 

/d d £ A 
(2Tr) d D„(k 1 )...D n (ki)D N (pi)...D N (pj) 

Here fi is a scale factor and A is the numerator function, which may contain anything. D w (k) = (£+ k) 2 — + ie 
and Dpj(p) = (£ + p) 2 — m 2 + ie are respectively the pion and nucleon propagator denominators, mo^ and m are the 
unrenormalized pion and nucleon masses appearing in the original Lagrangian. Thus the number of subscripts it and 
N correspond to the number of pion and nucleon propagators respectively. We will always redefine the integration 
variable I so as to make the first pion momentum k\ — and will drop it from the argument list. 

In general we can always reduce A to factors which can be removed from the integral or to powers of which at 
the one loop level can be reduced out using standard tensor expansions. Thus the only integrals which need to be 
evaluated explicitly have A = l. 

For the basic calculation we need the integrals ^[1], I-mr[q, 1], In\p, 1], lNN[Pi,Pf, 1], Ittn[p, 1] 5 I^NN[Pi,Pf,M, 
q,pi, 1]. For those diagrams with a mass insertion on an internal nucleon diagram, which duplicates one 
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of the nucleon propagators, we require the additional integrals Inn[p,P, 1], lNNN[Pi,Pf,Pf, 1], Innn [Pi,Pi,Pf, 1], 
I-kN n [P: Pi 1], I^NNN[Pi,Pf,Pf, 1], lTT-KNN[-<l,Pi,Pi, 1], where p can be either p, 4 or p 7 and where q=Pf~Pi- 

The evaluation of these integrals in the form needed for the IR or EOMS schemes proceeds in the standard fashion, 
as described for example in 12]. The meson and nucleon propagators are separately combined using the Feynman 
parameter approach. The two pieces are then combined and the infinities extracted using standard dimensional 
regularization formalism. The results can then be expressed in d dimensions in terms of the R and e = (4 — d)/2 of 
Eq. (|56[> and of relatively simple integrals over the Feynman parameters. This approach however leads, as discussed 
in the main text, to results which do not obey the usual HBChPT counting rules. 

Becher and Leutwyler ^| modify this procedure by dividing the integrals into two parts, one containing the infrared 
singularities and the other a regular polynomial in the expansion parameter. The regular part is then 'dropped', i.e. 
in a formal sense absorbed in the LEC counterterms. 

In order to discuss both the standard and the Becher-Leutwyler approach simultaneously we define a parameter 
j3 which flags the terms to be dropped in the Becher-Leutwyler procedure. Thus integrals involving only nucleons 
obtain an overall j3 IR in accord with the result that they are regular. Integrals with only pions are evaluated in the 
standard approach and so contain no j3 IR . For those integrals involving both pions and nucleon propagators the basic 
integral on the parameter z is divided into two parts and evaluated in accord with [Til IT2I ] as 

/■oo />oo 

dz -> J dz - I3 IR J dz . (A2) 

As discussed above, the singular terms proportional to R, which appear in both regular and infrared parts, can be 
recombined (i.e we eventually put j3 IR R — > R) and the renormalization of the LEC's carried out in the usual way. 
Thus (3 IR will serve to flag the regular terms which would be dropped in the Becher-Leutwyler procedure. 

With these preliminaries recorded we can list the results for the integrals we need in this calculation. 

/tt[1], Imr[q^ 1], In[p, 1], and InnIpiiP/, 1] are standard and the results are given here for completeness only: 

For/^1] 

m = ^k(R + H^))- (A3) 



For 1^ [q, 1] 



with 



(4tt) 2V v fi 



1 ,„ . . , ,mi 



W<?, 1] = (R + 1 + M^f )) + W„Aq 2 } > (A4) 



1 r 1 

Wnn[q 2 ] - ^2 J o dzMD-ir,) 



2 4 

q q 



9&Tr 2 mQ n 9607r 2 mQ W 



(A5) 
(A6) 



where here 



FovI N [p, 1] 



For I NN \pi,pf, 1} 



Q 2 

D = 1 - -V«(l - *) . (A7) 



" 



? 9 

m , , ,m 



I N [ P , 1] = I N [0, 1] = p IR j^(R + ln(^r)) • (A8) 



2 



I NN \p hP f,l] = l3 IR -^{R+l + ln(^j)) + W NN [q z }, (A9) 
I NN \p,p,l] = WO, 0,1] = f3 IR j^(R+l + ln(^)), (A10) 

lNNN\Pi,Pf,Pf,l] = lNNN\Pi,PhPf, 1] = W NNN [q 2 ] , (All) 
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where here 



with 

i r 1 

W NN [q 2 } = P IR -r^ / dzMp-iri) (A12) 
( 47r ) Jo 

= (j ^96ir 2 m 2 N + 96ir 2 m% ~ 9607r 2 m^ + ""^ (Al3) 

W NNN [q 2 } = / / (A14) 

(47r) 2 m 2 J a D - ir] 

327r*m^ 62ir z mj f 192TT 2 m]^ 

D = l-^z{l-z) and £ m = m " ~ m * , (A16) 
m 2 

For simplicity we have given only the first few terms in the expansions of the W's above, as the full expressions used 
are quite lengthy. In the actual calculations we kept more terms, as many as necessary to obtain the final amplitude 
through the first four orders in the expansion parameter. 

The remaining integrals involve both pion and nucleon propagators. For those we follow and generalize the procedure 
used in [12( for I^n- 

For In 7v [p, 1] we find, after combining denominators and evaluating via dimensional regularization, 

I wN \p,l] = -J—(R-l+ln(^))(2e-l) f dz . . (A17) 

Here C = C + C\{z - z ) 2 . For this case C\ = l + 2aQ + a 2 , C = a 2 (l- fi 2 )/Ci and z = a(a + Q)/Ci. Here (and 
below) a = m^jm and £1 = (p 2 — m 2 — mg T )/(2mmo T ). These integrals depend on the square of the four momentum 
p 2 but we will need them only at the physical on shell point p 2 = m 2 N . We must account for the fact that m ^ mjv 
and hence as above use m 2 — m 2 N = b m rri^ v where 8 m is a dimensionless parameter presumedly of order one. In fact, 
from Appendix iBl we have S m — 8ci?7ijv + This allows us to expand in powers of about the physical mass 
rriN- The integral can be done analytically, basically by integrating by parts, as in [l2T | and we obtain 



IMP, 1] = 7Tv» OS™ " ttt(P IR !)(! - U) + W„ N \p\m\mU , (A18) 
(4tt)^ 2m% 



where on shell 

\Al » r \mn 2 _ mn 2 mn _ 1 ( u , , — , , - _ _ 

167T 2 p 2 167TTOJV 327r 2 m 2 v 



ftl R 2 2 

W vN [m 2 N ,m 2 ,ml n ] = ^-M^) - 1) + -p*- + - W IR 



2 2 



+S m (p IR - 1) + (1 - S m )(\n(^) - /3 JK ln(^f ))) + (A19) 

In a similar fashion we find for I v NN\Pi,Pf, l],IwNNN\Pi,Pf,Pf, 1], and I w NNN\Pi,Pi,Pf, 1] 

1 m 2 f 1 C 1 z dz 

= -^ W («+ 1 + M^))(2e) ^ ^ ^ (c _ jT?)1+ « . (A20) 

InNNN\Pi,Pf,Pf, 1] = ^irJVJVJv[Pi,Pi,P/)l] 

1 f 1 Z" 1 z 2 dz 

= -sijspf + *>./„ •*/. P^F^ (A21) 

where now on shell with p 2 = p 2 = m 2 N , C\ = 1 + 2afi + a 2 — y(l — y)q 2 /m 2 , C = a 2 (l — fl 2 — y(l — y)q 2 /m 2 )/Ci 
and z — a(a + Q)/Ci. The integrals on z can be obtained analytically by generalizing the procedure of 0|. This 
result is then expanded in powers of the small parameter and integrated term by term on y. We thus obtain 

I W NN\Pi,Pf,l] = W^nnIp^p 2 ^ 2 ,™ 2 ,™^}, (A22) 

I*NNN\Pi,Pf,Pf,i\ = IffNNN[Pi,Pi,PfA] = ^—Wff NNN [p 2 ,p 2 f ,q 2 ,m 2 ,ml^}, (A23) 

m 07I 
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where on shell 



U — [nr{ ■ m{ ■ '"^ = + 3^r(^ln(^) - ln(^)) 



(i-U + ^r(i-2Ai-2g 



64^™^, m 647r 2 m^f 
64ir 2 m% 192w 2 m% [ ^to 2 / 

/3 ZR ln(^|)) + (A24) 

A* 



2 



• ! '"'^ (l-65 ro + --*+.... (A25) 



102471-774 ^ m ™> 512^^ 



The ttNN integral with duplicate nucleon propagators can be obtained by taking q 2 — > in InNN\pi-,Pf, 1], namely 

I 7r A'jv[p,P, 1] = W ff /viv[m^,m^,0, to 2 , m^] . (A26) 

Finally 



= --^{i+24 (A28) 



where we now have on shell with pf = to 2 ^, Cj = 1 + 2aO + a 2 — 3j(l — x)q 2 /m 2 , Co = a 2 — CiZq — x(l — x)q 2 /m 2 
and zo = a(a + Q)/C\ — a;(l — x)q 2 /(m 2 Ci). [44| Again we can do the z integration analytically and then expand in 
powers of the small parameter and do the x integration term by term. We then obtain 

I™N[-q,PiA] = — ^— W nnN \p 2 ,p 2 f ,q 2 ,m 2 ,ml^] , (A29) 

ImtNN[-1iPuPiA] = -^W^NNiphp^q 2 ,™ 2 ,™^], (A30) 

where on shell 

■,2 



W^ N [rn 2 Nl m 2 N ,q 2 im 2 ,m 2 ^] = — (1 + 9 ) 

327TTO/V 12777-0^ 



-^V(2 - «5 m - 3/?'* + ln(^L) - ^(^J) 

3271^777^ TO^ /J Z 

/,2 



-(1 + <W + , (A31) 



1927r 2 TO 2 v TOo 7 r 

W^m 2 mm l, q \m 2 ,rnU = -^ir(l + A- ) + (1 ^- > 



+ ^3 (! + 3* m ) + (A32) 

153D7rm^ r TOo„- 

Finally we observe that there is an alternative method for obtaining the loop integrals involving two nucleon 
propagators of the same momentum. It follows from the relations 

^ { ^ } { ^ } 2 ^ 19 (A33) 

g m 2 p2 _ m 2 p2 _ m 2 ' C?TO 2 TUq^ dS m 
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that one can get a loop integral with a duplicate propagator by taking derivatives, namely 

1 d 



m§ w d5 m 



1 d 

2 ax lNN\Pi,Pf,M = lNNN\Pi,Pf,Pf,M+ lNNN\Pi,Pi,Pf,M 

1 d 



ra 2 ^ dS m 
1 d 



lTTNN\pi,Pf, 1] = l7TNNN[Pi,Pf,Pf, 1] + lTrNNN[pi,Pi,Pf, 1] , 
2,r7riv[-?,Pi, 1] = /im-JViv[— 5,Pt,Pi,l] • 



(A34) 

(A35) 
(A36) 
(A37) 



We have checked that our results satisfy these relations. Note that to obtain I^nn \p, P, 1] for example to 0(m*) we 
need I v n[p, 1] to 0{m%) because of the m 2 introduced in the denominator by the derivative. This somewhat lessens 
the utility of this method for actually calculating the integrals with duplicate propagators. 



APPENDIX B: MASS AND WAVE FUNCTION RENORMALIZATION 

In the meson sector the pion mass and wave function renormalizations are calculated in standard fashion and are 
given in a number of sources, for example, In our conventions and notation we have 



m~ = m 0w 



1 ' ^( 



4(4tt) 



1 , / m 
In 



and 



Z* = l- 



2mi 
F 2 



1 lnf^ 



3(4tt) 2 3(4tt) 2 V M 



where the LEC's have been renormalized as 



R 



4(4tt) 2 
R 



(4tt) 2 



The renormalization of the pion decay constant is also standard and given from 5] by 

Ftx = Fa 



(4tt) 2 



(Bl) 

(B2) 

(B3) 
(B4) 

(B5) 



In the pion-nucleon sector the nucleon mass and wave function renormalizations must be calculated in a fashion 
consistent with the rest of the calculation. The appropriate diagrams contributing to the nucleon self energy are given 
in Fig. and the amplitudes corresponding to those figures are given by 

(B6) 
(B7) 
(B8) 

(B9) 
(BIO) 

+ m)^ 75 ]*. (Bll) 









M 2 NN 


= i*(4ciTO^)*, 




M?,NN 


= i*(2m^(8e 38 + 


ens + ene))* , 


M 4NN 


= -j^^iMpJ 


75 (^+ f+m) , 




Sim 2 ,- — . 
= ^*(2 Cl - 




M 6 nn 


3*ff4 / , 2 
= A p2 ( 4C 1™0. 


WUnn[p,P, ¥+ m)(jt 



We evaluate these amplitudes using standard dimensional regularization and expand the integrals in powers of the 
small momentum, keeping for now the finite terms which do not obey counting. As discussed above those finite terms 
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which would be dropped in the Becher-Leutwyler procedure are flagged with the symbol (3 IR . Likewise we flag terms 
from each diagram which would be dropped in the EOMS procedure with (3 EOMS . 

A first rcnormalization of the low energy constants is required to ensure that the difference between the physical 
nucleon mass mjv and the nucleon mass in the chiral limit mn is finite. In particular we take 

c\ = ci + lt^ 1 - 12 ^)' (B12) 

3i? 

e Ii5 + e ne + 8e 38 = eii5 + eii 6 + 8e 38 + ^^-^(c 2 - 8ci +4c 3 - 168cimjv.9i) • (B13) 
The finite renormalizations required read 

C IR = c EOMS ^ (B14) 

q 2 2 

IR , IR , rs IR _ EOMS , EOMS , o EOMS ggl m N9A f n i 1 ^W-^ 2 /„ 2 ^ 

e ii5 + e ii6 + 8e 38 - e U5 +e 116 + 8e 38 — — j-g— (28 + Lbin(m N /fj, )) 



3.9 



128F27r 2 mAr 



(3 + lnK/M 2 )), (B15) 



and 



cfOMS = c r _ 3mN9 A {1 + 8cimN _ (1 _ ^dm^Jin^/A* 2 )) , (B16) 

EOMS i EOMS i EOMS _ r ,r , or 3c 1 2 m jv g 2 1 oi„c„,2 /,,2\\ m 17 l 

e 115 + e 116 + e 38 - e 115 + e 116 + 8e 38 QF2--2 I 8 + «> m ("W^ jj . (B17J 

Ox , 7r 7r 

This leads to 



--H$(^ + M- 2 /m 2 )) 



-4™M K -2mt(e(f 5 + e(f 6 + 8 e ^) 



32F> 64F%n 2 m N K v 128F 2 tt 2 
"" r r ( 3c2 -32 Cl + 12 C3 )ln(m2/ M 2 ) + ^^^( M ). (B 18) 



64^2 V" ' L < 7T / / ^2 

Observe that all of the terms in the above formula which are proportional to 1/F% do not contribute to any physical 
amplitude we calculate here. They would lead to 1/F% corrections to loop diagrams, the same order as two loop 
contributions which we are neglecting. Furthermore the tree level diagrams do not contain nucleon masses for which 
these terms are relevant. Thus for practical purposes we can take 

m —> m N + Aciml + 2mJ(en 5 + en 6 + 8e 38 ) , (B19) 

where we have dropped all terms with 1/F%, which also allows {c{ R , e{ R 5 , e[ R e , e^ R } — > {ci, en5, en6, e 38 }. 
Note also that in the chiral limit only the R term survives, i.e. we obtain 

™°^= m - + M™ 2 //* 2 )) . (B20) 

It is common to introduce explicitly or implicitly a counterterm in the Lagrangian to eliminate the correction term 
in this formula and thus to interpret m from the beginning as the nucleon mass in the chiral limit, mjv, rather than 
as a 'bare' mass as we have done. 

Finally the nucleon wave function renormalization becomes 

Z " = l-^^(2(l+(3 IR ) + (3-2(3 IR )R + 16c imN (3 IR ((3 EOMS -l) 

+ 31n(4) - 2/J J *(l - 6c imN (l - f3 EOMS ))\K{™f) - ^) (B21) 
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If we reexpress the formulas for ttin and Zjq in terms of m and mo T and take (j, = m and take /r = then these 
results agree with those of Ref. [1^. Note however that there are no LEC's or counterterms available to absorb the 
(3 IR or p EOMS terms in Zjq, i.e. those terms which are to be dropped in the IR or EOMS schemes. Instead what 
happens is that these terms enter the amplitudes via the \fZ~N terms which appear in Eqs. I|15|) . i|16[l . (|17l) . (|45l) and are 
absorbed in LEC's elsewhere in the calculation. 
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FIG. 1: Diagrams which contribute to the coupling of external vector and axial vector currents to the nucleon. The solid, 
dashed, and wiggly lines correspond respectively to nucleons, pions, and external vector or axial vector fields. The unlabeled 
vertices come from C^ N whereas the ones labeled 2,3,4 come respectively from C^, C^ N and C^ N . 




\ / \ / \ / 
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FIG. 2: Diagrams which are higher order and which would not be included in an explicit calculation. However parts of these 
diagrams would be included implicitly by using the physical mass mjv in the propagators of loop integrals. 
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FIG. 3: Diagrams which contribute to the 7rNN vertex. 
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FIG. 4: The general loop integral 
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FIG. 5: Diagrams contributing to the nucleon self energy 



